Geometry and the Hidden Order of Luttinger Liquids: the Universality of Squeezed 

Space. 
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We present the case that Luttinger liquids are characterized by a form of hidden order which is 
similar, but distinct in some crucial regards, to the hidden order characterizing spin-1 Heisenberg 
chains. We construct a string correlator for the Luttinger liquid which is similar to the string 
correlator constructed by den Nijs and Rommelse for the spin chain. We re-analyze the spin one 
chain, introducing a precise formulation of the geometrical principle behind the so-called 'squeezed 
space' construction, to demonstrate that the physics at long wavelength can be reformulated in 
terms of a Z2 gauge theory. Peculiarly, the normal spin chain lives at infinite gauge coupling 
where it is characterized by deconfinement. We identify the microscopic conditions required for 
confinement thereby identifying a novel phase of the spin-chain. We demonstrate that the Luttinger 
liquid can be approached in the same general framework. The difference from the spin chain is 
that the gauge sector is critical in the sense that the Luttinger liquid is at the phase boundary 
where the Z2 local symmetry emerges. In addition, the 'matter' (spin) sector is also critical. We 
evaluate the string correlator analytically for the strongly coupled Hubbard model and further we 
demonstrate that the squeezed space structure is still present even in the non-interacting fermion 
gas. This adds new insights to the meaning of bosonization. These structures are hard-wired in the 
mathematical structure of bosonization and this becomes obvious by considering string correlators. 
Numerical results are presented for the string correlator using a non-abelian version of the density 
matrix renormalization group algorithm, confirming in detail the expectations following from the 
theory. We conclude with some observations regarding the generalization of bosonization to higher 
dimensions. 

PACS numbers: 64.60.-i, 71.27. +a, 74.72.-h, 75.10.-b 



I. INTRODUCTION 



The Luttinger liquid, the metallic state of one dimen- 
sional electron matter, is an old subject which is believed 
to be fully understood. In the 1970's the bosonization 
theory was developed which has a similar status as the 
Fermi-liquid theory, making possible to compute long 
wavelength properties in detail with only a small number 
of input parameters^. In the present era, the theory is 
taken for granted, and it has found many applications, 
most recently in the context of nanophysics 3 . Here we 
will attempt to persuade the reader that there is still 
something to be learned about the fundamentals of the 
Luttinger liquid. 

In first instance, it is intended as a clarification of some 
features of the Luttinger liquid which appear as rather 
mysterious in the textbook treatments. We make the case 
that a physical conception is hidden in the mathematics 
of the standard treatise. This physical conception might 
alternatively be called 'hidden order', 'critical gauge de- 
confinement' or 'fluctuating bipartite geometry'. It all 
refers to the same entity, viewed from different angles. 
This connection was first explored in our previous letter—, 
here we expand on these ideas to yield some practical 
consequences: (a) we identify symmetry principles al- 
lowing a sharp distinction between Luttinger liquids and 
for instance the bosonic liquids found in spin-1 chains 13 
(the 'no-confinement' principle, sections ITU IIII|) . (b) we 
identify a new competitor of the Luttinger liquid (the 



manifestly gauge invariant superconductor, section lTTTI a 
close sibling of the superfluid t — J model of Batista and 
Ortia&), and (c) these insights go hand in hand with spe- 
cial 'string' (or 'topological') correlation functions which 
makes possible unprecedented precision tests of the ana- 
lytical theory by computer simulations, offering also ad- 
vantages for the numerical determination of exponents 
(section EU). 

This pursuit was born out from a state of confu- 
sion we found ourselves in some time ago, caused by 
a view on the Luttinger liquid from an unusual angle. 
Our interest was primarily in what is now called 'stripe 
fractionalizationlZi2i2iiS. Stripes refer to textures found 
in doped Mott-insulators in higher dimensions. These 
can be alternatively called 'charged domain walls—: the 
excess charges condense on d — 1 dimensional manifolds, 
being domain walls in the colinear antiferromagnet found 
in the Mott-insulating domains separating the stripes. 
Evidence accumulated that such a stripe phase might be 
in close competition with the high- Tc superconducting 
state of the cuprate o 7 ' 12 and this triggered a theoreti- 
cal effort aimed at an understanding of stripe quantum 
liquids. The idea emerged that in principle a supercon- 
ductor could exist characterized by quantum-delocalized 
stripes which are however still forming intact domain 
walls in the spin system. Using very similar arguments 
as found in sections ITTI and ITTT1 of this paper, it can then 
be argued that several new phases of matter exist gov- 
erned by Ising gauge theory. This is not the subject 
of this paper and we refer the interested reader to the 
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Hterature£i2i2iiS. However, we realized early on that these 
ideas do have an intriguing relationship with one dimen- 
sional physics. 

Specifically, we were intrigued by two results which, 
although well known, do not seamlessly fit into the 
Luttingcr liquid mainstream: (a) the hidden order in 
Haldane spin chains as discovered by den Nijs and 
Rommelse^, (b) the squeezed space construction as de- 
duced by Woynarovichi^, Ogata and Shiba— from the 
U — ► oo Bethe Ansatz solution of the Hubbard model. 
As we will discuss in much more detail, after some fur- 
ther thought one discovers that both refer to precisely 
the same underlying structure. This structure can be 
viewed from different sides. Ogata and ShibaM empha- 
size the geometrical side: it can be literally viewed as a 
dynamically generated 'fluctuating geometry', although 
one of a very simple kind. Den Nijs and Rommelse ap- 
proached it using the language of order—: a correlation 
function can be devised approaching a constant value at 
infinity, signaling symmetry breaking. The analogy with 
stripe fractionalization makes it clear that it can also be 
characterized as a deconfinement phenomenon in the lan- 
guage of gauge theory. 

Whatever way one calls it, this refers to a highly or- 
ganized, dynamically generated entity. The reason we 
got confused is that there is no mention whatsoever in 
the core literature of the Luttinger liquid of how these 
squeezed spaces etcetera fit in the standard bosonization 
lore. To shed some light on these matters we found inspi- 
ration in the combined insights of den Nijs- Rommelse and 
Ogata-Shiba and we constructed a den Nijs type 'string' 
correlator but now aimed at the detection of the squeezed 
space of Ogata and Shiba. This is the principal device 
that we use, and it has the form, 



o str (\i-i\) = (Si n^(-i) 



Si 



(i) 



where Si is the spin-operator on site i while n; measures 
the charge density. By studying the behavior of this cor- 
relator one can unambiguously establish the presence of 
squeezed-space like structures. We spend roughly the 
first half of this paper explaining how this works and 
what it all means. In section [H] we start with a short 
review of the den Nijs-Rommelse work on the S = 1 
'Haldane— spin chains. This is an ideal setting to de- 
velop the conceptual framework. We subsequently re- 
formulate the spin chain 'string' correlator in a geo- 
metrical setting which makes the relationship with the 
Woynarovich-Ogata-Shiba squeezed space manifest. We 
finish this section with the argument why it is Ising gauge 
theory in disguise. This is helpful, because the gauge the- 
ory sheds light on the limitations of the squeezed space: 
we present a recipe of how to destroy the squeezed space 
structure of the spin chain. 

In section IIIII we revisit Woynarovich-Ogata-Shiba. 
The string correlator Eq. Q is formulated and subse- 
quently investigated in the large U limit. This analysis 
shows that the Luttinger liquid (at least for large U) 



can be viewed as the critical version of the Haldane spin 
chain. It resides at the phase transition where the gauge 
invariancc emerges, while the matter fields are critical as 
well. In this section we also argue why the squeezed space 
of the electron liquid cannot be destroyed. This turns out 
to be an unexpected consequence of Fermi-Dirac statis- 
tics. 

In the remaining two sections the string correlator 
is used to interrogate the Luttinger liquid regarding 
squeezed space away from the strong coupling limit. In 
section lTVl we demonstrate in a few lines a most surprising 
result: squeezed space exists even in the non-interacting 
spinful fcrmion gas! This confirms in a dramatic way that 
squeezed space is deeply rooted in fermion statistics; it 
is a complexity price one has to pay when one wants to 
represent fermion dynamics in one dimension in terms of 
bosonic variables. 

In section [V] we turn to bosonization. Viewing the 
bosonization formalism from the perspective developed in 
the previous sections it becomes clear that the squeezed 
space structure is automatically wired into the structure 
of the theory. In this regard, the structure of bosoniza- 
tion closely parallels the exact derivations presented in 
section lllll In section I VII we present numerical density 
matrix renormalization group (DMRG) calculations for 
the string correlators starting from the Hubbard model 
at arbitrary fillings and interaction strength, employing 
a non-abelian algorithm. These results confirm in a great 
detail the expectations build up in the previous sections: 
the strongly interacting limit and the non-interacting gas 
are smoothly connected and in the scaling limit the string 
correlator Eq. isolates the spin only dynamics regard- 
less the microscopic conditions. This has also practical 
consequences; we deliver the proof of principle that the 
non-universal exponents associated with the logarithmic 
corrections showing up in the spin-correlations can be 
addressed away from half-filling. From the combination 
of bosonization and the exact results for strong coupling 
suggest that the two point spin correlator can always be 
written in the scaling limit as the product of EqQJ and 
a charge-like string correlator, 



(S(x)S(O)) ~ D nn {x) O str (x) 



where the 'charge' string operator is defined as, 



(2) 



D nn (\i-j\) ee (n.(i) [lLj =i (-l) n " W ] n s (j)) (3) 

where n s (i) is 1 for a singly occupied site and otherwise. 
We confirm numerically that except for a non-universal 
amplitude the relation Eq. @ seems always satisfied at 
long distances. 

The conclusion to this paper addresses the broader per- 
spective including the relation to stripe fractionalization 
in 2+1 dimensions. 



3 



II. GEOMETRY, GAUGE THEORY AND 
HALDANE SPIN CHAINS. 

The 'Haldane-^ S — 1 Heisenberg spin chains are an 
ideal stage to introduce the notions of 'hidden order', 
squeezed space, and the relation with Ising gauge the- 
ory. These systems are purely bosonic, i.e. dualization 
is not required for the identification of the bosonic fields, 
and the powers of bosonic field theory can be utilized 
with great success to enumerate the physics completely. 
We refer in particular to the mapping by den Nijs and 
Rommelse^ on surface statistical physics. We are under 
the impression that this way of thinking is not widely dis- 
seminated and we start out reviewing some highlights. 
In the surface language, the meaning of 'hidden-order' 
becomes particularly simple (section III A|) . We subse- 
quently use these simple insights to reformulate this hid- 
den order in the geometrical language, the 'squeezed 
space' ( section III B|) . The next benefit of the Haldane 
chain is that the identification of squeezed space geome- 
try with Ising gauge theory is literal fsection lll C|) . This 
sets the conceptual framework within which we view the 
Luttinger liquid. 

A. Haldane spin chains: a short review 

Let us first review some established wisdoms con- 
cerning the Haldane spin chains. The relevant model 



is a standard Heisenberg model for S — 1 extended 
by biquadratic exchange interactions and single-ion 
anisotropy, 



H=Y / S l -S J +aY / (S>-S J ) 2 + DY / (SD 2 ■ (4) 

<ij> <ij> i 



In the proximity of the Heisenberg point (a = D = 0) 
the ground state is a singlet, separated by a finite en- 
ergy gap from propagating triplet excitations. It was 
originally believed that the long distance physics was de- 
scribed by an 0(3) non-linear sigma modeli^, suggesting 
that the ground state is featureless singlet. However, Af- 
fleck et aii& discovered that for a — 1/3, D = the exact 
ground state wavefunction can be deduced, having a par- 
ticularly simple form. This 'AKLT' wavefunction can be 
parametrized as follows. Split the S = 1 microscopic sin- 
glets into two Schwinger bosons 1-5 = 1, M s ) ~ b\ x a b\ 2a . 
The individual Schwinger bosons carry S — 1/2 and the 
wavefunction is constructed by pairing, say, the '1' bo- 
son with a '1' boson on the left neighboring site forming a 
singlet of valence bond, and the same with the '2' boson 
with its counterpart on the right neighbor, 



AKLT 



2 -N/2 



(tf frt it it i fit it 

l°i-l;lt°i;ll u i-l;ll u i-ip V y i;2T°i+l;2| 



X ( & !+1;1T 6 !+2;11 - & !+l;li 6 I+2;lT) ' ' ] l V&C ) 

I 



(5) 



This wave function clearly has to do with a transla- 
tion symmetry breaking involving nearest-neighbor sin- 
glet pairs, although in terms of spin-degrees of freedom 
which are different from the elementary spins. It has 
become a habit to call it 'valence bond solid order', 
i.e. to link it exclusively to the tendency in the spin 
system to form spin 1/2 singlet pairs. Den Nijs and 
Rommelse^ added a deep understanding of the physics 
of these bosonic spin chains by introducing the mapping 
on surface statistical physics. Although the AKLT wave 
function is a correct prototype for the ground state of the 
Heisenberg chain, it is not helpful with regard to what 
else can happen. On the other hand, by employing the 
formidable powers of surface statistical physics there are 
no secrets and it yields a natural view on the physics 
of the spin-chains. A highlight is their demonstration 
that this vacuum can be understood by a a non-local 
('topological') order parameter structure in terms of the 
real S = 1 spin degrees of freedom. The measure of 



order is the asymptotic constancy of a correlation func- 
tion. The conventional two-point spin correlator in the 
Haldane chain decays exponentially for large \i — j\, 



-\i-m 



(0) 



However, considering the following non-local spin cor- 
relator (or 'string' correlator), 



(st nU-i) s ' 



s z , 



constant 



(7) 



signaling a form of long range order which only becomes 
visible when probed through the non-local correlator Eq. 
0. For this reason it was called 'hidden order'. A main 
purpose of this section will be to introduce a more precise 
definition of this order. 

Den Nijs and Rommelse^ deduced the string correla- 
tor using the insights following from the path-integral 
mapping onto surface statistical physics. A first, crucial 
observation is that the natural basis for the spin chain is 



not in terms of generalized coherent states, but instead 
simply in terms of the microscopic M s = 0,±1 states 
of individual spins. Marshall signs are absent and these 
states can be parametrized in terms of flavored bosons 



(a) 



t ?,t 



u Q , u±1 subjected to a local constraint J2 a 



bLb. 



spin operators become, 



1. The 



S* = n i:1 - Hi -i 



(8) 



A second crucial observation is that because of the con- 
straint the problem is isomorphic to that of a (directed) 
quantum string living on a square lattice. This is some- 
what implicit in the original formulation by den Nijs and 
Rommelse, but used to great effect by Eskes et alw&. 
The mapping is elementary. A lattice string corresponds 
with a connected trajectory of 'particles' living on a lat- 
tice and this string can in turn be parametrized by a 
center of mass coordinate and the set of links connect- 
ing all particles. Consider only 'forward moving' links 
(the string is directed) and identify a nearest-neighbor 
link with a M s = bond, and an 'upward' and 'down- 
ward' next-nearest-neighbor link ('kinks') with M s = 1 
and M s = — 1 states of the spin on a site of the Haldane 
chain, respectively. It is easy to convince oneself that 
every state in the Hilbert space of the spin chain corre- 
sponds with a particular string configuration. The XY 
terms are responsible for the creation of kink-antikink 
pairs and the propagation of individual kinks along the 
string, while Ising terms govern the interactions between 
the kinks. In the path integral formulation, quantum 
strings spread out into world-sheets and the world-sheet 
of the lattice string corresponds with a surface statistical 
physics which is completely understood: the restricted 
solid-on-solid (RSOS) surface. 

The topological order of the Haldane spin chain trans- 
lates into a simple form of order in the surface language: 
the disordered flat phase. The +1 and —1 kinks on the 
time slice turn into up and down steps on the world-sheet 
in space time, see Fig. (QJ. In the disordered flat phase 
these steps have proliferated (kinks occur at finite density 
while they are delocalized) but on this surface every 'up' 
step is followed by a down step and the surface as a whole 
is still flat, pinned to the lattice. In the string language 
the order is therefore manifest, but it becomes elusive 
when translated back to the spin system. It implies that 
the ground state of the spin chain is a coherent superpo- 
sition of a special class of states. These are composed of 
indeterminate mix of 0, ±1 states. Take the 0's as a ref- 
erence vacuum and view the M s = ±1 states as particles 
carrying an in internal 'flavor' ±1. These particles are 
delocalized. However, every +1 particle is followed by a 
— 1 particle, modulo local violations (virtual excitations) 
which can be integrated out perturbatively (fig. [IJ. The 
hidden order is thereby nothing else than the staggered 
order of the ±1 flavors of the 'spin particles'. This or- 



_0_ _0 _0_ 

(b) 



FIG. 1: Mapping of the spin chain on a directed quantum 
string living on a lattice—. The M s = 0, ±1 states of the spin 
chain are equivalent to horizontal-, and upward/downward 
diagonal links tracing out the trajectory of the string on the 
lattice. The XY terms in the spin Hamiltonian correspond 
with the kinetic energy of the string problem causing both 
the creation of kink-antikink pairs (the ±1 bonds) and the 
coherent propagation of individual kinks. At the Heisenberg- 
and AKLT points hidden order is present. Although kinks 
are proliferated their direction is ordered: every up kink is 
followed by a down kink (a) . In the string representation this 
just means that the string pins to the links of the lattice. 
In the rough (XY) phase kinks have proliferated while their 
direction is disordered as well (b). 



der is not seen by the spin-1 operators because these also 
pick up the positional disorder of the ±1 'particles'. 



B. Squeezed space: sublattice parity as a gauge 
freedom. 

String correlators of the kind Eq. Q have the pur- 
pose of 'dividing out' the positional disorder with the 
effect that the order of the 'internal' ±1 flavors becomes 
observable. In order to see the similarity with the phe- 
nomena occurring in the Luttinger liquid we need a more 
precise description of how this 'division' is accomplished 
than that found in the original literature. It amounts 
to a geometrical mapping of a simple kind. The string 
correlator can be written in terms of the bosons as, 



(Sf Ui =i (-l) s i S?) = ((n 







x n^-i) 1 -"*." 



x (rij.i 



(9) 



n j.-i)) 



Why is this tending to a constant while the two-point spin 
correlator is decaying exponentially? From the discussion 
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in section III Al it follows that modulo local fluctuations the ground state wave function has the form, 
I 



|*) = ^ a(xx,x 2 , ...,x 2i ,x 2i+1 , ...) \x 1 (l),x 2 (-l), ...,x 2l (-l),x 2i+ i(l), ...) 

I 



(10) 



where the XiS refer to the positions of the ±1 particles 
on the chain, and the amplitudes a are independent from 
the 'internal' (±1) degrees of freedom; these 'internal' 
Ising spins show the antiferromagnetic order. In order to 
construct a two point correlator capable of probing this 
'internal' order it is necessary to redefine the space in 
which the internal degrees of freedom live. Start out with 
the full spin chain and for each configuration, whenever 
a site occupied by a 0-particle is found remove this site 
and shift, say, all right neighbors to the left, see fig. J2J- 
This new space is called 'squeezed space' and the effect of 
the map from 'full' to squeezed space is such that every 
configuration appearing in Eq. I|1U|) maps on the same 
antiferromagnetic order as realized on the squeezed chain. 



' n 
• 

B 
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p = +l 



= +1 
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B A 



p = +l 

FIG. 2: The geometrical mapping from 'full' (a) to 'squeezed' 
space (b). Given that some antiferromagnet lives in squeezed 
space, all that matters is the fate of the sublattice parity p. 
When sites are reinserted, the sublattice parity of the sys- 
tem in squeezed space flips every time a hole is passed when 
viewed in full space. These sublattice parity flips hide the 
order present in squeezed space. 

Obviously, if it would be possible to probe squeezed 
space directly, the hidden order would be measurable us- 
ing conventional two point correlators. The string cor- 
relator achieves just this purpose. All that matters is 
that the order in squeezed space is a staggered (antifer- 
romagnetic) order. For such order one needs a bipar- 
tite geometry: it should be possible to divide the lattice 
into A and B sublattices such that every site on the A 
sublattice is neighbored by B sublattice sites and vice 
versa. In one dimension any space is bipartite (even the 
continuum). This subdivision can be done in two ways: 
■■■A-B-A-B-- or ■ ■ ■ B - A - B - A- ■ corre- 
sponding with the Z 2 valued quantity we call sublattice 
parity. For a normal lattice the choice of sublattice parity 
is arbitrary, it is a 'pure gauge'. However, in the map- 



ping of squeezed to full space it becomes 'alive', actually 
in a way which is in close correspondence to the work- 
ings of a dynamical Z 2 gauge field as will become clear 
later. Consider what happens when squeezed space is un- 
squeezed (Fig|5}. When a '0' particle including its site is 
reinserted, the 'flavor' site, say, on its right side is shifted 
one lattice constant to the right. The effect is that rela- 
tive to the reference sublattice parity of squeezed space 
the sublattice parity in unsqueezed space changes sign 
every time a '0' particle is passed. The effect is that flips 
in the sublattice appear to be 'bound' to the particles 
viewing matters in full space. In order to interrogate the 
'flavor' order in squeezed space one has to remove these 
sublattice parity flips. This can be achieved by multi- 
plying the spin with a minus one every time a particle 
is encountered: (—1) x (— 1)' +1 = (-1)'. The den Nijs 
string operator is constructed to precisely achieve this 
purpose, 



(Ki - rn,-i) |n£ =i (-l) 1 -'H.°J (rij-,1 - 
= ((-ly^x-ni,^) [nf =i (-l)"».o 



nj-i)) 



A B 



A B 



B 



A B 



(11) 

Hence the string correlator measures the spin order in 
squeezed space by removing the sublattice parity flips. 
The positional disorder of the particles is equivalent to 
motions of the sublattice parity flips, scrambling the or- 
der living in squeezed space, and these are removed by 
the string operators. 

• • o • • • o • o • • 

ABABABABABA 



A 
1- 



FIG. 3: Same as Fig. but now for the case that an Ising 
long range order is present in squeezed space, corresponding 
with the hidden order of the S = 1 Heisenberg spin chain. 
From such pictures one can directly deduce the workings of 
string operators by just focusing on the spins. However, mat- 
ters are equally meaningful in the absence of long range order 
in squeeze space, and the formulation terms of a geometrical 
mapping involving sublattice parity is more general. 

The above argument emphasizes the geometrical na- 
ture of the mechanism hiding the order. It might at this 
point appear as a detour because one arrives at the same 
conclusion by just focusing on the 'flavor' orientations, 
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see Fig. @. However, as will become clear in later sec- 
tions, the construction is still applicable even when the 
spin system in squeezed space is disordered. Hence, it is 
more general and rigorous to invoke the geometrical sub- 
lattice parity as a separate degree of freedom in addition 
to the degrees of freedom populating squeezed space. 



C. Squeezed spaces and Ising gauge theory. 

At first sight, it might appear that sublattice parity 
is not quite like a normal dynamical degree of freedom. 
However, it is easily seen that it is nothing else than an 
uncommon ultraviolet regularization of Z2 gauge fields. 
From the above discussion it is clear that the 'flavor' 
degrees of freedom of the ±1 particles can be regarded 
as independent from their positions in unsqueezed space. 
These flavors are Z 2 valued and can be measured by, 



(1 



&lAo)(-i)'sf 



(12) 



The positions of the particles drive the uncertainty in the 
value of the sublattice parity and these are captured by 
the Zi valued operators, 

of = (— 1)™ !, ° (13) 

and it follows that modulo a factor of order 1, 



((ni,i 



m.-i) nt.(-i) 1 -"'' 



3,1 



,-1) 



and in the presence of the hidden order, 



(14) 



(15) 



i.e., at distances large compared to £ the correlations 
between the r spins have disappeared but they re-emerge 



n? =i of 



is attached to every 



when the operator string 
spin. 

This suffices to precisely specify the governing symme- 
try principle: the long distance physics is governed by 
a Z 2 gauge field (the o's) minimally coupled to spin- 1/2 

matter (the r's). The strings ILi =i af simply correspond 

with the Wilson loop associated with the Z2 gauge fields 
rendering the matter correlation function gauge invari- 
ant. The two point correlator in the r's is violating gauge 
invariance and has therefore to disappear. This gauge 
invariance is emerging. It is not associated with the mi- 
croscopic spin Hamiltonian and it needs some distance £ 
before it can take control. Therefore, the gauge-violating 

( T i T j) i s non zer0 f° r I* — j\ < £• 

This is an interesting and deep connection: the inde- 

terminedness of the sublattice parity in full space is just 
the same as invariance under Z2 gauge transformations. 
One can view the squeezed space construction as an ultra- 
violet regularization of Z2 gauge theory, demonstrating 



a simple mechanism for the 'making' of gauge symme- 
try which is distinct from the usual mechanism invoking 
local constraints (e.g., references ITollT^ ) 

Is this yet another formal representation or does it re- 
veals new physical principle? As we will now argue, the 
latter is the case. Viewing it from the perspective of 
the gauge theory, it becomes immediately obvious that 
there is yet another possible phase of the spin chain: the 
confining phase of the gauge theory. To the best of our 
knowledge this phase has been overlooked because its ex- 
istence is not particularly obvious in the spin language. 

For a good tutorial in gauge theory we refer to Kogut's 
review^. Focusing on the most relevant operators, the 
Z2/Z2 theory can be written as, 



Z = 



s 



VrVcre 



(fxdr 



7~i &i j Tj T" 



/ v UplaqO' 
plaq 



(16) 



leaving the gauge volume implicit in the measure, r and 
a are Z 2 valued fields living respectively on the sites and 
the links of a (hypercubic) space-time lattice. The ac- 
tion of the gauge fields is governed by a plaquette action, 
i.e. the product of the fields encircling every plaquette, 
summed over all plaquettes. The gauge invariance cor- 
responds with the invariance of the action under the flip 
of the signs of all the o's departing from a site i, ac- 
companied by a simultaneous flip of the Tj. This gauge 
invariance implies that tj = 1 <-> — 1 and ( n Tj ) = 
while ( tj pro] Tj ) can be non-zero (with T a line 
of bonds on the lattice connecting i and j ; the Wilson 
loop). This is the most general ramification of the gauge 
symmetry and Eq. (|15f) is directly recognized. 

The relation between the gauge theory and the 
squeezed space construction is simple (Fig. Q). The 
gauge invariance is just associated with the indetermi- 
nacy of the sublattice parity in unsqueezed space. If the 
O's would not fluctuate one could ascribe a definite value 
to the sublattice parity everywhere, and this is equivalent 
to choosing a unitary gauge fix in the gauge theory. How- 
ever, because of the derealization of the O's one cannot 
say if the sublattice parity is +1 or —1 and this corre- 
sponds with the gauge invariance. 

As is obvious from the string correlator, the Z2 gauge 
fields (coding for the indeterminacy of the sublattice par- 
ity) are coupled to matter degrees of freedom being just 
the 'flavors' living in squeezed space. In the hidden- 
order/disordered flat phase these are Ising spins showing 
long range order. The constancy of the string correla- 
tor at long distances reflects this fact. From the view- 
point of the gauge theory this appears as an absurdity. It 
means that the hidden order phase is the Higgs phase of 
the Z2/Z2 gauge theory, characterized by a gauged mat- 
ter propagator becoming asymptotically constant. In the 
gauge theory this can only happen in the singular limit 
where the gauge coupling K — > 00! 
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FIG. 4: Squeezed space mappings as geometrical interpreta- 
tion of Ising gauge theory. Although the word-lines (in space- 
time x,t) of the ±1 particles span up a bipartite lattice for 
an observer which is just watching word-lines, this bipartite- 
ness is hidden in full space when the particles are delocalized. 
This is equivalent to the conventional lattice regularization 
of a Zi gauge theory involving a plaquette action where the 
+ — gauge invariance of the link variables acquires the 
meaning that it is impossible to determine the bipartiteness 
of squeezed space by measuring in full space. The absence of 
free visons (minus fluxes) does imply that the hidden bipar- 
titeness exists and the existence of squeezed space corresponds 
with deconfinement. Taking the unitary gauge is equivalent 
to squeezing space. 



Even under the most optimal circumstances (high di- 
mensionality) , a Wilson loop should decay exponentially 
with a perimeter law due to local fluctuations in the 
gauge sector. Stronger, it is elementary that in 1+1 D 
the Higgs and deconfining phases are fundamentally un- 
stable to confinement. This law can only be violated in 
the singular limit K — > oo. Hence, the hidden order ap- 
pears as highly unnatural within the framework of the 
gauge theory. What is the reason that confinement is 
avoided in the Haldane spin chain? More interestingly, 
what has to be done to recover the natural confinement 
state? 

The disorder operators in the gauge sector are the vi- 
sons or gauge fluxes. These are point-like entities (in- 
stantons) in 1+1 dimensional space-time. For any finite 
value of the coupling constant K these will be present at 
a finite density with the result that the vacuum is confin- 
ing and the implication that ( n [n r er] Tj ) — > at large 
distances. Translating this to the geometrical language, a 
vison corresponds with a process where a squeezed space 
of even length on time slice r turns into a squeezed space 
of odd length on time slice r + 5t or vice versa. In this 
way a minus gauge-flux is accumulated on a time- like pla- 
quette (see Fig. |SJ). In terms of the degrees of freedom 
of the spin chain this means that a single M s = state 
can fluctuate into a M s = ±1 state and vice versa. It 
is obvious now why the spin chain corresponds with the 
K — > oo limit of the gauge theory, namely the Hamiltio- 
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FIG. 5: As Fig. (3) except that now a single ±1 particle is 
annihilated. This means that the bipartiteness of squeezed 
space is destroyed and this is in one-to-one correspondence 
with the presence of isolated visons (minus fluxes) in the Z2 
gauge theory formulation: the 'natural' confining state. 



nian of the former only contains pairs of spin raising or 
lowering operators ~ Sf SJ . From Eq. (JSJ) it follows im- 
mediately that '0' particles can only be created or anni- 
hilated in pairs. These processes do change the length of 
squeezed space but they turn even-length squeezed space 
into even-length squeezed space, or odd-length squeeze 
space into odd-length squeezed space. Confinement re- 
quires odd to even or even to odd fluctuations. In the 
geometrical language, deconfinement means that space- 
time is still bipartite although the two ways of subdivid- 
ing space-time are indistinguishable. Confinement means 
that bipartiteness is destroyed outright because squeezed 
space-time can no longer be divided in two sublattices 
due to the presence of the visons. 

Going back to the spin chain the remedy becomes ob- 
vious. First, one should impose a preferred direction of 
the spin quantization axis which can be accomplished 
by taking for instance D > 0. Apply now a transver- 
sal field: Sf = (B/2) £ 2 [S+ + fir] . Given Eq. 

(JSJ, it follows that for any finite value of B isolated O's 
will turn into ±l's and squeezed space can be destroyed! 
A somewhat delicate aspect is that the matter field is 
in the fundamental representation (i.e. the ZijZi mat- 
ter/gauge theory). We will analyze this elsewhere in more 
detail but it can be easily argued that for any finite trans- 
verse field a featureless 'Higgs-confinement' phase will 
take over, characterized by exponential decay of the den 
Nijs string correlator. 

One can actually wonder whether under experimental 
circumstances the hidden order phase can be ever truly 
realized in spin chains. The hidden order is in a sense 
pathological, as deconfinement is in 1+1D: a transverse 
magnetic field of any strength suffices to destroy the hid- 
den order. Is it at all possible to devise experiments 
such that transverse fields are rigorously vanishing? Ex- 
ploiting the relationship with gauge theory, a number of 
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other interesting conclusions can be reached regarding of 
the spin chains. However, spin chains are not the real 
subject of this paper, and we leave this for a future pub- 
lication. The primary aim of this section is to supply 
a conceptual framework for the discussion of the more 
convoluted 'hidden order' in the Luttinger liquids. Let 
us list the important lessons to be learned from the spin 
chains, and indicate how these relate to the Luttinger 
liquids: 

1 . The central construction is squeezed space, the ex- 
istence of which can be detected using den Nijs-type 
string correlators. The determination of such a cor- 
relator for the Luttinger liquid is the subject of the 
next section. 

2. The phases where sublattice parity flips are truly 
delocalized are characterized by an emergent Z2 
gauge symmetry. We make the case that such 
phases can in principle occur also in the Luttinger 
liquid context, while the Luttinger liquid itself re- 
sides right at the phase boundary where the 
local invariance emerges. 

3. In the spin chains squeezed space can be destroyed 
by transverse fields causing confinement. We argue 
that in the Luttinger liquids this is impossible be- 
cause of the fermion minus signs of the electrons, 
with the ramification that squeezed space is uni- 
versal. 



III. LUTTINGER LIQUIDS: SQUEEZED SPACE 
IN THE LARGE U LIMIT. 

The focus in this section is entirely on the Luttinger 
liquids which can be regarded as continuations of those 

I 



The charge part ipSF represents the wave function of non- 
interacting spinless fermions where the coordinates Xi re- 
fer to the positions of the N singly occupied sites. The 
spin wavefunction ipHeis. is identical to the wave function 
of a chain of Heisenberg spins interacting via an antifer- 
romagnetic nearest neighbor exchange, and the coordi- 
nates yj, j = 1, . . . ,M refer to the M positions occupied 
by the up spins in the Heisenberg chain. The surprise is 
that the coordinates yj do not refer to the original Hub- 
bard chain with length L, but instead to a new space: a 
chain of length N constructed from the sites at coordi- 
nates Xi, X2, xn given by the positions of the charges 
(singly occupied sites) in a configuration with amplitude 
ipSF- One immediately notices that it is identical to the 



describing the long distance physics of Hubbard models. 
The bottom line is that these Hubbard-Luttinger liquids 
are characterized by a critical form of the spin-chain type 
hidden order as discussed in the previous section. This 
criticality has two sides: (a) the (Z2) gauge fields are 
critical, in the sense that the Luttinger liquid is associ- 
ated with the phase transition where the local symmetry 
emerges, (b) the matter fields (spins) are also in a critical 
phase. 



The argument rests again on the squeezed space con- 
struction, and this should not come as a surprise to 
the reader who is familiar with the one dimensional 
literature. This construction was introduced first by 
Ogata and Shiba— , who rediscovered earlier work by 
Woynarovichli regarding a far-reaching simplification in 
the Lieb and Wu Bethe-Ansatz solution of the Hubbard 
model— in the U — > 00 limit. This Woynarovich-Ogata- 
Shiba work just amounts to the realization that in the 
large U limit the structure of the Bethe-Ansatz solution 
coincides with a squeezed space construction. For sim- 
plicity, assume a thermodynamical potential n > such 
that no doubly occupied sites occur. For U tending to 
infinity, the ground state wavefunction ip of a Hubbard 
chain of length L occupied by N electrons (with N < L) 
factorizes into a simple product of spin- and charge wave- 
functions, 



(17) 

I 

squeezed space mapping for the Haldane spin chains dis- 
cussed in the previous section, associating the M s = 
states of the spin chain with the holes and the M s = ±1 
states with the singly occupied sites carrying electron 
spin up (+) or down (— ). In fact, as already pointed out 
by Batista and Ortia&, one can interpret the spin chain 
as just a bosonic t — J z mode, i.e. lowering the SU (2) 
symmetry of the Hubbard model to Ising, dismissing the 
Jordan- Wigner strings making up the difference between 
spinless fermions and hard-core bosons, and last but not 
least adding an external Josephson field forcing the holes 
[M s = 0, in the spin language) to condense giving a true 
Bose condensate. 

Since the geometrical mapping is the same, a 'string' 



«/>(£!, ...,X N ;yi,.. .,Vm) = ^Sf{x\, ■ ■ .,X N ) ^iTeis.^l, ■ ■ ■ ,Vm) ■ 
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operator equivalent to that of den Nijs and Rommelse 
can be constructed for the Luttinger liquid. In order to 
measure the spin correlations in squeezed space starting 
from unsqueezed space one should construct an opera- 
tor which removes the sublattice parity flips. Define the 
staggered magnetization in unsqueezed space as, 



M(x) = (-l) x S(x) 



(18) 



Compared to the corresponding quantity in squeezed 
space, these acquire an additional fluctuation due to the 
motions of the sublattice parity flips. Since these flips 
are attached to the holes, they can be 'multiplied out' by 
attaching a 'charge-string', 



(M') z (x) = M z {x){ 



(i-«t°t(j)) 



(19) 



where 1 — ntot(j) is the number of holes on site j and the 
charge operator n tot {j) = n^{j) +n±(j) taking the values 
0, 1 and 2 for an empty-, singly- and doubly occupied 
site, respectively. (M') z is representative for the 'true' 
staggered magnetization living in squeezed space. The 
action of the charge string Hj(—l) 1 ~ ntot ^ is to add a 
— 1 staggering factor only when the site j is singly occu- 
pied, thereby reconstructing the bipartiteness in squeezed 
space. It follows that the analogue of the den Nijs topo- 
logical operator becomes, 

O str (x) = ( {M'Y{x) (M')«(0) > 

= (M z (x) {-l)^=° {1 ~ ntotU)) M z {0)} (20) 
= -(S z (x)(-l)^=i nt ° tU) S z (Q)} . 

The focus of the remainder of the paper is on the anal- 
ysis of this correlator. To the best of our knowledge, 
correlators of this form have only been considered before 
in the context of stripe fluids in 2+lD& 2 £. String correla- 
tors have been constructed before in the one dimensional 
contexljSiS^ but these are of a different nature, devised 
to detect 'hidden order' of an entirely different type. 

On this level of generality it might appear that the 
hidden order of the Haldane chain duplicates that of the 
Luttinger liquid. However, dynamics matters and in this 
regard the Luttinger liquid is quite different. Instead of 
genuine disorder in the 'charge' sector and the true long 
range order in the 'spin' sector of the spin chain, both 
charge and spin are critical in the Luttinger liquid and 
this makes matters more delicate. 

We learned in the previous section that in order to 
measure the hidden order one should compare the con- 
ventional two point spin correlator (M(r)M(O)) with the 
string correlator defined in Eq. (|20|l . Let us compute 
these correlators explicitly in the large U limit. In the 
calculation, the string correlator turns out to be a sim- 
plified version of the two point correlator. The latter 
was already computed by Parola and Sorella™ starting 
from the squeezed space perspective. Let us retrace their 
derivation to find out where the simplifications occur. 

Start with the observation that a Heisenberg spin an- 
tiferromagnet is realized in squeezed space. This implies 



that the squeezed space spin-spin correlator has the well- 
known asymptotic form, 



Heis .(j) EE {S*(j)S'(0)) 

_> (_iyr lnl/2(i) 



(21) 

(-iyo stag u) 

where L is a constant 26 , while j labels the sites in 
squeezed space. 

The charge dynamics are governed by an effective sys- 
tem of non- interacting spinless fermions. Define their 
number operators as n(l) where / refers to sites in full 
space. Define the following correlation function, to be 
evaluated relative to the spinless fermion ground state, 



PsfU) = (n(0)n(x)6 Vn(l) 



)SF 



(22) 



By definition this measures the probability of finding j 
spinless fermions in the interval [0,x], given one fermion 
located at site and one at site x. Parola and Sorella 22 
show that the exact relation between Eq. (|21|) and the 
two point correlator in full space is, 

x+1 

(S z (x)S z (0)) = J2 P SFU)°Heis.(j-l) 



x+1 



J=2 
x + 1 



-J2 p sFU)(-iyo stag u) 



Let us now consider instead the string correlator, 
O str (x) = -<S*( )(-l)£J=i" (j V(x)) 

,(j)(-i) J '- 2 CW(j 



(23) 



1 



(24) 



x + 1 

- J2 P SFU)Osta 9 U) ■ 

The difference between the two point correlator and the 
string correlator looks at first sight to be rather unre- 
markable. The staggering factor (— 1) J associated with 
the sign of staggered spin in squeezed space (Eq. |2"T1) sur- 
vives for the two point correlator, but it is canceled for 
the topological correlator because (— 1) J ~ 2 x (— l)-? -1 = 
(— l) 2j ~ 3 = — 1. However, this factor is quite important 
because it is picked up by the charge sector due to the 5 
function appearing in the definition of Psf (Eq. I22fl . 

In Eqn.'s 1)231241) spin and charge are still 'entangled' 
due to the common dependence on j. However, it can 
be demonstrated that asymptotically this sum factorizes. 
It can be proven 2 ^ that the sum SJ=2 P sf(J){~^-)' 1 f{j) 
with f(j) bounded and satisfying 



] - 3 



< 2P 



In 1 / 2 ! (a 



(25) 
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differs from the sum 

x+l 



E p ^(i)(-i) J 

J=2 



/«r>«) 



(26) 



where (pt 



N/L is the fermion density) 

1 -r+nw) 



^x+1 
(n(0)n(x))sF Z-i] 
x Ptot + 1 — > Xptot 



(27) 

by terms vanishing faster than log 2 ^ . The condition 
Eq. (|2*5|l is satisfied by the squeezed space staggered 
magnetization f(j) ~ O s t a g(j) and since the above re- 
sult does not depend on the presence of the staggering 
factor (— iy it applies equally well to the two point spin 
correlator and the string correlator. 

Given this factorization property, let us first consider 
the string correlator, 



(28) 



x+l 

O str (x) = Y. P SF{j)Osta g {j) 
J'=2 
x+l 

J^PsfU) O stag {x P tot) 
It is easy to demonstrate that the sum over the Psf is 



just the density-density correlator of the non-interacting 
spinless fermion system 



E-^Phti) = (n(Q)n(x)) SF 



Ptot 



cos{2kpx) 



(29) 



with hp — nptot- We arrive at the simple exact result, 



O str (x) = (n(x)n(0)} 



SF-J-ln^ip^ + O 
= Y^\^\p t0t x)+0{^1). 



ln 3 / 2 (x) 
72 



(30) 

This confirms the intuition based on the squeezed space 
picture. The topological correlator just measures the spin 
correlations in squeezed space which are identical to those 
of a Heisenberg spin chain, Eq. (|21|l . At short distances 
this is not quite true, but it becomes precise at large 
distances due to the asymptotic factorization property 
Eq. Of course, O str measures in units of length of 

the full space and because in squeezed space sites have 
been removed the unit of length is uniformly dilated x — > 
Ptotx. By the same token, the amplitude factor reflects 
the fact that there are only p to t spins per site present in 
full space. 

The calculation of the two point spin correlator is less 
easy. Using again the factorization property, 



(S'(x)S'{0)) = 



x+l 

x+l 

E p sfO")(-i) j a 

j=2 



ta 9 {{r)x)+0(^L 



(31) 



-Dnn,SF\X) T h O \ — 



PtotX 



r 



Due to the staggering factor, the 'charge function' 
Dnn,SF{x) is now more interesting, 



x+l 



J=2 
x+l 



^2{n(0)n(x)S[J2n(l)-j))sF 



i=2 



vi=0 



= („(0) (-l)£r=o n[l) n{x)) 



SF 



(32) 

The spin correlations are modulated by a function re- 
flecting the uncertainty in the number of sublattice parity 
flips which can be expressed in terms of expectation val- 



ues of charge string operators. For spinless fermions the 
following exact identity holds for the number operator, 



n(j) = \ [l - (-!)"« 



(33) 



which implies, 

D nn , S F = \ [D SF {x - 2) + D SF {x) - 2D SF (x - 1)] 

(34) 

demonstrating that this function is the second lattice 
derivative of the charge-string correlator, 



D SF (x) = <(-l)5Xo n «) 



SF ■ 



(35) 
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Even for free spinless fermions this function has not 
been derived in closed analytic form. However, it can be 
easily evaluated numerically and we show in the appendix 
that it is very accurately approximated by, 

( (_l)£J=i "Ci) \ sp = ^ 2 %/2 cos(7rp tot (x- 1)) 
^sin(irp tot ) \/x - 1 

(36) 

where A is a constant evaluated to be A — 
0.6450002448 22 . Using Eq. O it follows immediately 
that, 

D nn , SF (x) = (n(x) (-l)K nU} n(0)) SF 
A 2 (cos(itpsf) — 1) cos(ttpsfx) 

y/2sm(np SF ) \[x ( 37 ) 

A 2 cos(2kpx) 

where, as before, 2kp = irptot and introducing the charge 
stiffness K c which takes the value 1/2 in a free spinless 
fermion system. This is the desired result, and combining 
it with Eq. (|31|l we arrive at the asymptotically exact 
result for the two point spin correlator in the large U 
limit, 

{S'(x)S'(0)) = A 2 V2T COS{ l k + F K X) l^ 2 (x/2) 

O X c 

This calculation demonstrates quite explicitly why the 
spin correlations in this Luttinger liquid are sensitive to 
the charge fluctuations. The latter enter via the uncer- 
tainty in the location of the sublattice parity flips which 
is expressed via the function D nn or equally the more 
fundamental function D. Due to the factorization prop- 
erty Eq. (|31|l it enters in a multiplicative fashion. The 
string correlator is constructed to be insensitive to the 
sublattice parity fluctuation and it follows that, 

{S*(x)S'(0)) ex ^(^(x)(-l)K nt ° t(j) ^(0)). 
x n ° 

(39) 

This is in close analogy with Eq. I|15l) for the Haldane 
chain. The difference is that in the spin chain the string 
correlator is decaying exponentially slower than the two 
point correlator while in the large U Luttinger liquid the 
difference is only algebraic. This has an obvious reason. 
In the spin chain, the 'charge' sector is truly disordered 
(Bose condensed), such that the 'charge-charge' corre- 
lations decay exponentially and this will obviously also 
cause an exponential decay of the charge string correlator 
D. The charge sector in the Luttinger liquid is critical, 
exhibiting algebraic correlations. As have demonstrated 
explicitly above, this also renders D to be algebraic. We 
argued in section[H]that the exponential difference found 
in the spin chain signals the emergence of an Ising gauge 
symmetry: the charge string just corresponds with the 



Wilson loop of the gauge theory. By the same token, the 
algebraic difference in the Luttinger liquid means that 
the Ising gauge symmetry is not quite realized. How- 
ever, power laws indicate criticality and this is in turn 
associated with a second order phase transition. Thus 
we are considering a correlator which measures directly 
the gauge fields; its power law characteristic indicates 
that the gauge symmetry itself is involved, and the logi- 
cal consequence is: the Luttinger liquid is located at the 
continuous phase transition where local Ising symmetry 
emerges. 

This sounds odd at first sight. However, one should 
realize that this Ising gauge symmetry is just dual to the 
superfluid phase order in the charge sector. Although in 
1+1D true long range superfluid order cannot exist, the 
Luttinger liquid can be viewed as an entity which is at the 
same time an algebraic superfluid and an algebraic charge 
density wave. In principle, when one applies an external 
Josephson field acting on the charge sector alone it will 
directly turn into a true superfluid. In this superfluid 
the number correlations are short ranged and this implies 
that the charge-string will decay exponentially. 

A caveat is that this Josephson field has to be applied 
in such a way that the spin system is unaffected. For 
instance, applying a standard Josephson field acting say 
on the singlet channel ~ VPf 1 !^ has the automatic effect 
that a spin gap opens and one can continue adiabatically 
to the strong singlet pairing limit. At long distances, 
all that remains is doubly occupied sites and holes and 
it is no longer possible to construct squeezed space. It 
is 'eaten' by the spin gap. However, at least in prin- 
ciple one can construct a 'charge only' Josephson field. 
Consider the large U limit. The Bethe-ansatz wavefunc- 
tion demonstrates that the ground state in the decoupled 
charge sector is in one-to-one correspondence to that of 
a free spinless fermion Hamiltonian. One can simply add 
to this Hamiltonian a Josephson field acting directly on 

the spinless fermions ~ Hj J2<ij> c \ c j an d for any finite 
strength of Hj the charge ground state will correspond 
with a BCS superconducting state. By construction, this 
field will leave the squeezed space structure and the spin 
sector unaffected. The ramification is that the quantiza- 
tion of number density is truly destroyed and since holes 
continue to be bound to the sublattice parity, the disor- 
der in the number sector becomes the same as Z2 gauge 
degeneracy in the spin sector. This is the same type of 
construction as suggested by Batista and Ortia& in their 
identification of the Haldane spin chain with a superfluid 
t — J z model. 



IV. SQUEEZED SPACE AND 
NON-INTERACTING ELECTRONS. 

The existence of squeezed space is remarkable, and in- 
tuitively one might think that one needs highly intricate 
dynamics associated with strong electron-electron inter- 
actions in order for squeezed space to have a chance to 
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emerge. The evidence for its existence presented so far is 
entirely based on very special strongly interacting cases 
(the Haldane spin chain, the large U Hubbard model) 
which can be solved exactly for more or less acciden- 
tal reasons. However, in the previous paragraphs we 
have constructed and tested a measuring device which 
can unambiguously detect squeezed space also in cases 
where simple exact wave functions are not available. Al- 
ternatively it can be detected even in cases where one 
knows the wavefunction but where the squeezed space 
structure is deeply buried because the coordinates are 
not of the right kind. Our measuring recipe is straight- 
forward: compute the string correlator Eq.J3U|) and find 
out if behaves like the pure spin chain, or whatever 'mat- 
ter' system one expects to populate squeeze space. 
The simplest possible example is the non-interacting, 



spinful electron system. As we will demonstrate using 
only a few lines of algebra, it survives the test! We in- 
terpret this as a remarkable feat of the fermion minus 
signs. Squeezed space refers eventually to a bosonic rep- 
resentation of the fermion problem, and apparently the 
minus sign structure in terms of the fermion representa- 
tion is of sufficient complexity to make possible an entity 
as organized as squeezed space in the boson language. 

The proof is as follows. For a system of S = 1/2 
fermions we can use the following operator relations, 

S z (y) = U n t(y)- n l(y)) (40) 
ntot{y) = ni(y)+ni(y) . 

The string correlator can be written as, 



O str {x) = - (S*(x) (-l)^=i nt ° t{j) S*(0)) 

= -i<„ T(a;) (-i)£;:No)^ 

[411 



In the non- interacting limit, the spin up and spin down 
electrons behave as two independent species of free spin- 
less fermions. Since the expectation value of any operator 



involving only either up- or down spin creation and an- 
nihilation operators is the same, Eq. 141|) simplifies to, 



O str = -(S z (x) nt ° tU) S z (0)) = 

_i <7lSF(x)( _ 1) E;: 1 1 »-o')^ (0 ))((-i)E; = ->-(i) ) (42) 



where the operators now refer to spinless fermions. We 
recognize in this expression the D$ F and the D nn g F 
we already encountered in section ITTT1 (Ean. 's I32I35|I . In 
addition we also need, 

D n , S F = (n SF (x)(-l)^ nsF{j) ) (43) 
= | (D SF (x - 2) - D SF (x - 1)) 

being the first lattice derivative of D, using once again 
the operator identity Eq. ll.'j.'SII . The topological corre- 
lator can therefore be expressed entirely in terms of the 
'fundamental' string operator Ds F (x) ~ (n(— 1)" SF ) as, 

O str {x) = \ [D SF (x - 2)D SF (x) - D SF (x - l) 2 ] (44) 

o 

The function Ds F (x) was already encountered (Eq. 



1361 see also the appendix) and using this result, 
O str (x) = illHtal 

_ A 4 sin(fc F ) (45) 

A 4 sin(7rp to t/2) 
— Ax ' 

where p$ F = ptot/2 — [p] +Pi)/2- Note that 2k F = 
nptot = TT 2N y F and so k F = 7T ^ F = irp SF . We also 
calculated the string correlator numerically using the 
method explained in the appendix. 

In figureElwe show the numerical result for O s tr{x) for 
a density p to t = 2N SF /V = 0.2 and V = 200 which is 
in excellent agreement with the analytic expression Eq. 
I)45[). In Figure we show the numerical results for vari- 
ous densities on a log-log plot highlighting the algebraic 
decay with an exponent K s = 1. 

It is obvious where this exponent, equal to unity, is 
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FIG. 6: The function O s tr = -{S'(x)(-l)^>i=i ntotU) S z (0)) 
for (7 = calculated numerically using the algorithm dis- 
cussed in appendix A. Here p tot — 2Nsf/V = 0.2 and 
V = 200. The drawn line is the analytic solution Eq. 14511 . 
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FIG. 7: The function O str = -<S z (z)(-l)^=i ntotU) S* (0)) 
as computed numerically for the free spinful fermion gas at 
densities p tot = 0.2, p to t = 0.6 and p tot = 1, shown in a log- 
log plot. The algebraic decay implied by Eq. 14511 is indicated 
by the straight line. 



coming from in the calculation. From Eq. 144|) it follows 
that O s tr ~ l/(x C|SF ) 2 where the spinless fermion expo- 
nent K Cj sf = 1/2- This looks at first sight rather unspec- 
tacular but one has to realize that the two point spin cor- 
relator of the free-fermion gas decays faster, (SS) ~ l/x 2 
and the topological correlator therefore uncovers a more 
orderly behavior. Furthermore, the only symmetry rea- 
son to expect such an exponent to be equal to unity is 
the protection coming from SU (2) (spin) symmetry. Can 
we be certain that this result proves that even in the 
non-interacting limit a Heisenberg chain is lying within 
squeezed space? The above computation is not very ex- 
plicit in this regard and the persuasive evidence is still to 
come: bosonization, and especially the numerical results 
presented in section IVT1 showing that the asymptotic be- 
havior of the string correlator is independent of U and 



density. 



V. SQUEEZED SPACES AND BOSONIZATION. 

Arriving at this point, we are facing evidence that the 
squeezed space is actually not at all special to the large 
U limit. It could well be ubiquitous in one dimensional 
electron systems. How does bosonization fit in? After all, 
during the last thirty years overwhelming evidence accu- 
mulated for bosonization to be the correct theory in the 
scaling limit. Squeezed space is of course fundamental; it 
is among others a precise description of the meaning of 
spin-charge separation. How could bosonization ever be 
correct if it would not somehow incorporate the squeezed 
space structure? In section IV Bl we make the case that 
the peculiarities in the structure of the theory, originat- 
ing in the core of the bosonization 'mechanism' (i.e., the 
Mandelstam construction for the fermion operators), are 
just coding for squeezed space. Again, the string opera- 
tor is the working horse. By just tracking the fate of the 
string- and two point spin correlators in the bosonization 
framework, it becomes evident that it is in one-to-one 
correspondence with the strong coupling limit. This ob- 
servation is further amplified in Appendix B where we 
discuss an intuitive argument by Schulz which turns out 
to subtly misleading. To fix conventions, let us start out 
collecting some standard expressions. 



A. The bosonization dictionary. 

To fix conventions let us collect here the various stan- 
dard bosonization expressions we need latest. At the 
Tomonaga-Luttinger fixed point the dynamics is de- 
scribed in terms of gaussian scalar fields tp s and tp c for 
spin and charge, respectively. Introducing conjugate mo- 
menta Tl s . c the Hamiltonian is, 



TL 



fi=c,s 2 



Jdx 



K^l + Md^^f , (46) 



where K s (v s ) and K c (v c ) are the spin and charge stiff- 
ness (velocity) respectively. For globally SU (2) symmet- 
ric spin systems K s = 1 and K c is depending on the 
microscopy, but generally < K c < 1 for repulsive inter- 
actions. 

Electron operators can be re-expressed in terms of 
these bosonic fields via the Mandelstam construction. 
Starting from the spinful Dirac Hamiltonian describing 
the linearized electron-kinetic energy, 

(47) 

The field operators of the left- (Vv) and right (ipcr) 
moving fermions are expressed in terms of the bose fields 
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as, 



1p<r(x) 



i^/Tr(<p(x)-f dyU(y)) 



(48) 



where r\ a , fj a are the Klein factors keeping track of the 
fermion anti-commutation relations. 

Starting from the normal ordered charge density the 
total charge density can be written as, 



n to t{x) = : njjx) + nj. (x) 

2 Oipc 



^ + C dw(x) + OI dw (x), 



(49) 



where d x ip c represents uniform components of the charge 
density, while the various finite momentum components 
are lumped together into Ocdw- The dominant contri- 
butions come from momenta q — 2kp and 4fcp, 

Ocdw(x) = 2 k F (x) + 4kp (x) 



e -2«F* e n/2^c(z) cos [^2^(2:)] 



O ikF {x) 



(50) 



Similarly, the spin operator S z {x) becomes, 

— :nf(x)-n l (x): 

= \f^9 x (p s (x) + S dw,z(x) + 0\ DWz {x) , 

(51) 

where d x ip s refers to the uniform (ferromagnetic) compo- 
nent while the finite wavevectors are dominated by the 
q = 2ftF component, 

°SDW,z( x ) - ®Sz,2k F (x) 



2tt 



In addition we need the usual rules for constructing the 
propagators of (vertex) operators in a free field theory 
like Eq. (pEHj). 



(d x (p ll (x)d x (p v (0)) 

/ g in V27T [ipu, (x)- ipp (0)] ^ 



r 27r ^ (53) 

3 



-A",, 



B. Vertex operators and squeezed space. 

It is a peculiarity of bosonization that the charge field 
enters the spin sector in the form of a vertex operator 



~ e^ c , see Ea. l|5"2"|) . This can be traced back to the 
Mandelstam construction for the fermion field operators, 
Eq. I)48[l. indicating that the fermions are dual to the 
fields ip: the fermions have to do with solitons or kinks 
in the bose fields. 

Let us observe the workings of bosonization from the 
viewpoint offered by the strong coupling limit discussed 
in section m We found that the charge-string correlator 
D(x) is the most fundamental quantity keeping track of 
the fluctuations in the sublattice parity. Let us see what 
bosonization has to say about this correlator. 

This function becomes in the continuum, 



D{ X ) = ((-i)E;= 

^Y, X j=0 n tot{j) 

7r J*dy n tot (y) 



{cos 
(cos 



(54) 



The theory is constructed to represent the scaling limit 
and therefore we should focus on the leading singulari- 
ties. According to En.q49[l. the total charge is given by 
n to t — x/^Jft d x (p c plus finite q components. One can 
easily convince oneself that the latter will give rise to 
subdominant contributions which can be neglected in the 
scaling limit. Hence, 



D(x) 



(cos 



(cos 
l 



7T f*dy {y^pK d y cp(y) 



(55) 



Since bosonization can only probe non-zero wave vec- 
tor components of the density the expressions are correct 
up to multiplicative factors ~ cos(wpx) (p is average den- 
sity). Keeping this in mind, the outcome is fully consis- 
tent with the result obtained for the large U case (Eq. 
H3 K c = 1/2 in this limit) but now extended to arbitrary 
values of the charge stiffness! 

The correspondence between bosonization and the 
strong coupling analysis becomes very obvious in the 
derivations of the two point spin correlator and the 
string correlator. Let us recall the standard derivation 
in bosonization of the spin correlator, 
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{S>(x)S*(0)) 
{OsdwAx)0 ] sdw,M) 



1 ,d(p s (x) dip s (0), 



2ir dx 



dx 



(O 



SDW,z\ 



(x)Ol DW M) + h.c. 



1 

8^ 
1 

8^' 



e -2ik F x ^ e iv / 27r[^c(a:)- ¥ >c(0)]^ e 'iv / 27r[¥> s (2;)-V s (0)]\ 



(56) 



-2ikpx 



l 

r K c +K s 



and the spin-spin correlation function becomes 



(S'(x)S'{0)) 



K s 1 



1 cos{2kpx) 

4^ 



r K c +K s 



Comparing this with the large U outcome, Eq.'s 1|H8[I. 
the correspondence is clear: ^\^Wb{x)-^ b (0)\\ > j g ^ e 
staggered magnetization of the spin chain in squeezed 
space, Eq. J5lJ. In strong coupling, the sublattice par- 



ity fluctuations enter via the function (n(-l)2> n ) (Eq. 
I37fl which differs from D by just a factor cos(7rpx). This 



(57) is of course precisely e 



-2ik F x ^ e iV2Tr[^ c (x)~<p c (0)]^ j n 

bosonization expression Eq. (|56|l . Notice that the sub- 
dominant uniform component ~ 1 /x 2 was just ignored in 
the strong coupling analysis. 



The correspondence is further clarified by considering 
the string correlator. Straightforwardly, 



O str (x) = -(^(x)(-l)K ntot(j V(0)) 

= fi(4^(x)e*^(^^(°«9 rfs (0)) + /..c.) 

/ e-^F* ^ e iV2^(yc(x)-y e (0)) e -»^(y e (a)-y e (0))^ e -»V2^(y s (x)-y s (0))^ + h ( ,\ (58) 



And these contributions add up to, 



11 1 cos(2fe F a;) 1 cos(2fei?2;) ^ ' 

~~ ~~ 47T 2 x 2 + K " ~ 4-rr 2 x 7 ^ 4ir 2 i K >+ 4K t ' 



The first term is obviously the (over corrected) uni- 
form magnetization and the leading singularity at finite 
wavevectors is, 

(S z (z)(-l)K ntot0 V(0)> = C0S(2 ^ } . (60) 

Again the caveat applies that bosonization cannot keep 
track of the average charge density and the oscillatory 
factor in the numerator should therefore be ignored - 
this 'flaw' is just inherited from D(x), Eq. j55jl. Where 
is this leading singularity coming from? It corresponds 
with the third line in Eq. 1)58(1. This algebra is expressing 
that the charge vertex operator coming from the charge 
string exactly compensates for the charge vertex opera- 
tors attached to the spin operators. We recognize that 
this is in precise correspondence with Eq. 's (|24I30|) of the 
strong coupling limit. The charge string is coding for the 
fluctuating kinks in the sublattice parity and the string 
correlator is constructed to remove these from the spin 
correlations. 



What have we achieved? The above leaves no doubt 
that the algebraic structure of bosonization is exactly 
coding for the structure we discussed in a geometrical 
language in section ITTT1 However, in section ITTT1 we had 
to rely on the simplifications arising in the strong cou- 
pling limit. The algebraic structure of bosonization is 
however universal and independent of microscopic con- 
ditions like the strength of U. For instance, in the non- 
interacting limit K c — K s — 1 and one directly infers 
that the bosonization expressions Eq. (|rJUl ISTJl are con- 
sistent with the exact results we derived for the string- 
and spin correlators for this limit in section HVI Although 
there are some caveats regarding the use of bosonization 
to calculate (charge) string correlators, these are entirely 
of a technical nature and these affect only subdominant 
singularities: see appendix B. We can therefore safely 
conclude that bosonization is just encoding the squeezed 
space geometrical structure which is manifest in strong 
coupling. The 'hard-wired' structure of bosonization, in 
combination with the string operators, leaves no room for 
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any other conclusion that squeezed spaces are ubiquitous 
in Luttinger liquids. It is indeed the case that even non- 
interacting one dimensional electron systems have deep 
connections with hidden order in Heisenberg chains. 

VI. NUMERICAL RESULTS 

To verify that the correlator O s t r indeed demonstrates 
that squeezed space exists for finite values of the Hub- 
bard coupling U/t and arbitrary density, we performed 
numerical calculations using the DMRG method 2 ^. The 
DMRG is an ideal tool for these purposes, because the al- 
gorithm construction implies that string correlators are, 
in principle, no more difficult to construct than ordi- 
nary two-point correlators. Indeed, the string opera- 
tor (— l)" s is precisely that which is already used to en- 
sure the correct commutation relations for the creation 
and annihilation operators. We utilized the non-abelian 
formulation^ of the DMRG, which makes use of the 
SU{2) <g) SU{2) ~ 50(4)/Z 2 spin and pseudo-spin sym- 
metry of the Hubbard model™, thereby giving a substan- 
tial improvement in efficiency. The pseudo-spin symme- 
try is an expansion of U(l) particle number symmetry N 
to an SU (2) symmetry which we denote here by Q (this 
is sometimes also denoted by /). In the 50 (4) represen- 
tation, the particle-number is given by the z-component 
of the pseudo-spin, N — 2Q Z + 1. In our calculation, the 
basis states are 50(4) multiplets, labeled by two half- 
integral quantum numbers (s,q) denoting the total spin 
and total pseudo-spin respectively. 

Addressing the scaling limit with the DMRG method 
is subtle. In the DMRG method, the ground-state wave- 
function is calculated in a Hilbert space which is trun- 
cated. The parameter controlling the truncation is the 
number of states kept in each 'block', to. The actual 
dimension of the space in which the ground-state wave- 
function is determined is of order (4m) 2 . This trunca- 
tion introduces an error which, for a 'well-behaved' sys- 
tem, is completely systematic and can be corrected for 
by calculating the appropriate scaling as m — > oo. For 
the ground-state energy, this scaling is understood and 
a routine calculation in DMRG. For correlation func- 
tions, the scaling is highly non-linear and difficult to per- 
form, not least due to a result highlighted by Ostlund 
and Rommer— : the wavefunction obtained by DMRG 
is a (position-dependent) matrix-product wavefunction, 
which implies that the long-range asymptotic behav- 
ior of all two-point correlation functions is exponential, 
with a correlation length that depends on the number 
of states kept m. While in principle one can determine 
this correlation length and fit the remaining (algebraic) 
components of the correlation function, this is in fact 
not necessary due to a not so well understood property 
of (position-dependent) matrix product wavefunctions, 
namely in the short-distance correlations the exponential 
due to the finite truncation is not present at all. Thus, 
as long as a sufficiently large number of states are kept 



to be close to the scaling limit at distances less than the 
characteristic transition point where the correlator be- 
comes exponential, the exponents of algebraic terms can 
be determined with high accuracy without any additional 
corrections due to the finite truncation. 

Also of note is that matrix-product wavefunctions 
generically carry long-range string order, in the sense 
that it is likely that all string correlation functions decay 
exponentially in the asymptotic limit, but it is permissi- 
ble that the decay is to a non-zero constant. The canoni- 
cal example is the AKLT wavefunction, which is obtained 
exactly in (non-abelian) DMRG with m = 1 states kept. 
In principle, the variational nature of DMRG implies that 
for a finite number of states kept one could inadvertently 
and incorrectly obtain a state that has non-zero string or- 
der. This is not a serious issue and is entirely analogous 
to the case of ordinary two-point correlators which, in the 
absence of a symmetry constraint, may have a spurious 
(but usually negligible) non-decaying component. For ex- 
ample, a not-quite-zero uniform magnetization resulting 
in a non-zero constant in the spin-spin correlator. The 
point is that the construction of DMRG treats hidden 
order of the den Nijs-Rommclsc type on a very similar 
footing as more conventional order. 

In the calculations presented here, we used m — 1000 
50(4) states kept, and a lattice size of L = 1000. The 
lattice size was chosen to be rather large in an attempt to 
reduce the effect of the open boundary conditions. How- 
ever this is not strictly necessary and the usual averaging 
procedure suffices to eliminate the Friedel oscillations and 
obtain the correct scaling form of the correlators even for 
much smaller lattices. 

We calculated the string correlator O s t r , Eq. the 
sublattice parity correlator D, Eq. (|35|) and its second 
lattice-derivative D nn , Eq. ©, for a large variety of 
filling factors p = 0.1... 0.9 and U/t = . . . 16. No- 
tice that the number operators appearing in the 'charge' 
strings D and D nn correspond with n s measuring the 
presence (1) or absence (0) of a singly occupied site. In 
the exponent one might as well take the total charge 
density n tot — + rtj., i.e. (—1)™*°* = (-1)™=. How- 
ever, D nn ~ (n s Il(-l) n >n 8 ) + (n tot II(-l) n ">«n tot > be- 
cause n s cannot distinguish empty from doubly occupied 
sites whereas n tot does. On the bosonization level this 
subtlety does not matter, but it is consequential for the 
numerically exact charge string correlators. As the strong 
coupling analysis in section ITVl demonstrates, the charge 
string coding for the squeezed space structure is actu- 
ally D nn because empty and doubly occupied sites are 
indistinguishable in the squeezing operation. 

The obtained correlation function O s t r appears in Fig. 
©, plotted on a log- log scale. It is clear from the figure 
that the leading order term in O s t r is algebraic, with 
an exponent that is independent of both the filling factor 
and U . The fitted exponent is equal to 1, with a variation 
over all parameter ranges of ~ 5%. We percieve this as 
a striking result^, taking away all doubts regarding the 
'universality of squeezed space': regardless microscopic 
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circumstances we have identified a correlation function 
which always behaves as if the electron system is just the 
same spin-chain. 

Even the small variation of the exponent is explainable, 
employing logarithmic corrections. At the U/t — > oo 
Woynarovich-Ogata-Shiba point, the wavefunction fac- 
torizes exactly and the O s t r correlator measures exactly 
the logarithmic corrections of the isotropic 5 = 1/2 anti- 
ferromagnetic Heisenberg chainSHiSI. This coincides with 
the well-known form at half-filling2£, where the presence 
of the charge gap implies Heisenberg-like behavior of the 
logarithmic corrections for any U > 0. However, as far as 
we know, for finite U /t away from half-filling, the exact 
form of the logarithmic corrections is not known. It is 
well understood that these corrections are not governed 
by the conformal field theory, implying that these are 
non-universal quantities which are sensitive to the short 
distance dynamics. Hence, they are not necessarily in- 
dependent of doping and interaction strength. This is 
confirmed by the non-interacting limit U = where the 
system has no log corrections, and there is no reason 
to expect a discontinuity between U = and U — ► oo. 
Therefore we take the scaling form of the string correlator 
to be 



O str (x) = A{p, U) 



ln Q (x) 



(61) 



and fit for the running exponent a. The obtained ex- 
ponent appears in Fig. 0. We emphasize that this is 
a very rough calculation, obtained by a direct fit of the 
finite-size data to the asymptotic form, ignoring finite- 
size corrections. For the Heisenberg model these correc- 
tions are important 3 ? therefore our U — > oo estimate of 
a ~ 0.25 — 0.3, deviates from the exact value of a = 1/2. 
Indeed, our result is rather reminiscent of earlier Heisen- 
berg model calculations 34 which suffer from similar is- 
sues. A careful scaling analysis, done by Hallberg, Horsch 
and Martinez for the Heisenberg chainS, should present 
no difficulty and will be reported in a subsequent paper. 
However, from the present results we can already safely 
conclude that a is a function of U jt and perhaps also 
density. 

We have argued in previous sections that the charge 
fluctuations present in the ordinary two-point correlators 
are due to sublattice parity fluctuations. We found in 
section lHll that in the strong coupling limit the following 
rigorous result holds for the staggered component of the 
spin-spin correlator, 



(S(O)S(x)) ~ O str (x)D nn (x) 



(62) 



Our argument is that bosonization reflects this structure 
and we are now in the position to test this relation nu- 
merically for arbitrary values of U and density. As we 
already emphasized, to isolate the squeezed space the 
number operators in D nn should measure the density of 
singly occupied sites, n s . In addition, away from the 
Woynarovich-Ogata-Shiba point Eq. I|62|l is not longer 
exact but it should become exact in the scaling limit. 
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FIG. 8: The Ostr correlator on a log-log plot, for a wide 
variety of filling and couplings. The slope, which determines 
the exponent (up to log corrections) of the leading order term, 
is equal to 1 independent of the parameters. 
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FIG. 9: The parameter a governing the logarithmic correc- 
tions to Ostr- The logarithmic corrections are very sensitive 
to finite-size effects, which were not taken into account here. 
Thus the strong coupling limit deviates from the expected 
a = 1/2. 



Eq. I|62l) should hold up to a U, p dependent prefactor 
factor which is set by short distance physics. This is ex- 
actly what we find. This is demonstrated by Fig. 110|) 
which shows the exponent of the D nn correlator, which 
turns out to be given by 



D nn {x) = B(p, c/) C05( ^ FX) + 0(x-^) 



(63) 



where K c is the usual density and [/-dependent charge 
stiffness of the Hubbard model. It follows that 



(S(Q)S(x)) = F(p, U) C ° S }^ ] In- 1 , S 



(64) 



coincident with the well known asymptotic behavior of 
the two point spin correlator in the Luttinger liquid. This 
completes our case. The fact that we not only isolate the 
spin-only dynamics in the Luttinger liquid using O s t r but 
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that we can reconstruct the two point spin correlator by 
dressing it with an entity which is exclusively counting 
the sublattice parity mismatches (D nn ) leaves no doubt 
that squeezed space is universal. 

Let us end this section with giving some numerical 
results regarding the non universal prefactors A(p,U), 
B(p,U) and F(p,U). These are clearly sensitive to the 
details of the short wavelength dynamics and have there- 
fore a similar status as non-universal amplitudes in any 
critical theory. Hence, these have to be calculated nu- 
merically. 
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FIG. 10: The exponent of D nn (x). This function isolates the 
charge contribution to the correlation functions, hence gives 
a direct determination of K c . The solid lines are guides to 
the eye. 



The prefactor of the O s t r string correlator, A(p, U) is 
given in fieurc lTTl The numerical prefactor coincides with 
the expected exact expression at U — and follows the 
expected form cx p for U — > 00, for a Heisenberg chain 
diluted by a hole density of (1 — p). The exact slope of the 
U — > 00 prefactor depends sensitively on the exponent of 
the log corrections. The underestimation of a in equation 
(|61[l results in the prefactor of figure ITT1 being somewhat 
large; the U — > 00 form should be^ 



p\n 1/2 {px) 



(2tt) 3 / 2 



(65) 



This differs from the correlator of a stretched Heisen- 
berg chain by a prefactor p 2 , which is due to the dilution 
of the spins; for the Heisenberg chain (s) — 1/2, but for 
the Hubbard model (s) — n s /2. Thus, with all prefac- 
tors accounted for, the factorization of the U — * 00 spin 
correlator ia 22 i 35 



(S(0)S(x)) = --±— + ±O str (x)D nn (x) 

3 3A 2 cos(2fc F ) cos(2fc F a;)ln 1/2 (2;) ( 66 ) 

~ ~A(^xJ + (2tt)3/ 2 p^sin(2fc F ) = ' 



with 2fc F = Ttp. 

For finite coupling the exact factorization of the wave- 
function is destroyed by local fluctations, so Eq. l|rjrj|) 
only applies rigorously in the strong coupling limit. As 
shown in section IIVI however, the scaling form applies 
even to U — 0, with the introduction of a non-universal 
amplitude T(U,p), 

(S(0)S(x)) = --JL—+T(U,p)O str (x)D nn (x) . (67) 

Figure 1121 shows this amplitude as a function of den- 
sity and U, which is always finite implying that squeezed 
space is ubiquitous. 



VII. CONCLUSIONS: THE FERMION MINUS 
SIGNS. 

In first instance the pursuit presented above can be 
seen as an exploration of the usefulness of string corre- 



lators of the den Nijs and Rommclsc type in the con- 
text of one dimensional physics. To our perception these 
correlation functions are worthy additions to the stan- 
dard repertoire of one dimensional physics. This will be 
further amplified in a next paper where we will further 
explore the information one can obtain from string cor- 
relators like D and D nn . 

In this paper we used string correlators to clarify some 
conceptual issues in one dimensional physics. String cor- 
relators go hand in hand with the simple geometrical 
ideas which emerged in the study of Haldane spin chains 
and the strong coupling Bethe Ansatz solution of the 
Hubbard model. These correlators make it possible to 
address to what extent these notions are of relevance 
to generic Luttinger liquids and we made the case that 
squeezed spaces are hard-wired into Luttinger liquid the- 
ory. It is merely a matter of recognition. 

Although complementary to the standard descriptions, 
we find that the squeezed space notion does exert unify- 
ing influences. It is not an accident that we started out 
discussing the Haldane spin chains. We hope that we con- 
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FIG. 11: The prefactor of O s tr- The numerical data at U = 
matches the exact form determined in section ITVI The U — > 
oo prefactor is proportional to the density, exactly as required 
for a diluted Heisenberg chain. 
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FIG. 12: The aplitude F obtained from (S(0)S(x)) ~ 
YO s tr(x)D nn (x). For clarity plotted as Fp 2 , which is equal 
to unity in the strong coupling limit. The heavy solid and 
dashed lines are the exact expressions at U = oo and (7 = 
respectively. The light dashed lines linking the numerical data 
are guides for the eye. 



vinced the reader that there is a unity underneath which 
becomes obvious in this language, while it is far from 
obvious in the standard formulation of bosonization. 

Is it more than just clarification? If so, it should be 
that these insights can be used to deduce states of one di- 
mensional quantum matter which have been overlooked 
before. In the Luttinger liquid context we have deduced 
one such novel state: the 'charge only' superconductor we 
introduced at the end of section ITTT1 This entity can also 
be discussed in the bosonization language. It is a prereq- 
uisite to drive the system away from critically such that 
the charge sector is genuinely disordered. This requires 
an external Josephson field stabilizing superfluid phase 
order. A conventional Josephson field acting on electrons 
pairs in the singlet channel is expressed as (recall section 



Hj = Bj J dx (V' T (x)V'x(x) - ipi(x)ipf(x)) 
~ Bj J dx cos [\/27r#c] sin [\/27n/?s] 



(68) 



involving the dual charge field d x 9 c (x) — — n(x). This 
imposes phase order (pinning of C ) but it has also the 
immediate effect of opening a spin gap (~ sin [ ^/2Trip s ~j ) . 
This spin gap means that the spins are paired in pairwise 
singlets and a squeezed space cannot be defined for these 
singlets. Instead, what is required is a Josephson field 
acting exclusively on the charge fields, 



H' T = B.i 



dx cos 



2TT0 r 



(69) 



This will enforce disorder on the charge sector, leaving 
the spin sector unaffected. Recalling the discussion of the 
spin chain, this charge disorder turns into a Z 2 gauge in- 
variance in the spin sector. The spin system in squeezed 
space resides at the (SU(2)) critical point separating the 
XY and Ising fixed points and together with the minimal 
coupling to the deconfining Z 2 gauge fields a state of mat- 
ter is realized which is symmetry-wise indistinguishable 
from the critical state of the Haldane spin chain found at 
the transition from the hidden-order phase to the 5 = 1 
XY phase. 

Although such a state is a theoretical possibility, it is 
less clear whether it can be realized in nature. Bosoniza- 
tion is helpful in clarifying this issue. Starting out with 
electron operators, it appears to be impossible to con- 
struct a Josephson field of the form Eq. (|69|) . One will al- 
ways find that the charge Josephson field is accompanied 
by a (relevant) operator in the spin sector. This might 
well turn out to be a fundamental obstruction. In the 
one dimensional universe the charge and spin fields are 
more fundamental than electrons, and a-priori Eq. I|69|) 
is physical. However, a Josephson field will in practice 
correspond with a mean-field coming from three dimen- 
sional interactions and this implies that this mean-field 
has to be a composite of electron degrees of freedom. 

As we argued, squeezed space is hard-wired into the 
bosonization formalism and even exotic states like those 
discussed in the previous paragraphs are in principle 
within the reach of the formalism. By implication, if a 
state of electron matter would exist where squeezed space 
is destroyed, it would be beyond bosonization. In the 
context of the (bosonic) spin matter of the Haldane chain 
we encountered this possibility. Helped by the identifi- 
cation of the Z 2 gauge symmetry, we presented a recipe 
(the transversal field) to stabilize a non-squeezed space 
('confining') phase of the spin chain. Is this also possible 
in the electron Luttinger liquids? 

In this regard it is helpful to view these matters from 
a yet another angle: the Marshall signs introduced by 
Weng in the one dimensional context as an addition 
to the squeezed space construction needed to describe 
fermion propagators; see also reference for the ex- 
tension to 2D and for some interesting observations re- 
garding Marshall signs and spin-charge separation in ID. 
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Marshall signs refer to the theorem that the ground state 
wave function of a S = 1/2 spin system defined on a 
bipartite lattice with nearest neighbor exchange interac- 
tions is nodeless: it is a bosonic state. In the strong cou- 
pling limit the spin system in squeezed space is of this 
kind, and this explains in turn why the Bethe-Ansatz 
solution reveals that the charges are governed by spin- 
less fermions. The total wavefunction has to be anti- 
symmetric and because squeezed space exists the spin 
sector is symmetric, so that the fermionic grading resides 
in the charge sector. 

Although we are not aware of an explicit proof, it has to 
be that this 'division of statistics' is universal in the scal- 
ing limit. Our string correlator demonstrates that at long 
distances the squeezed space spin system does behave ex- 
actly like the (unfrustrated) Heisenberg chain and it is 
hard to imagine that this would survive a drastic change 
involving the nodal structure of the spin wavefunction. 
Let us assume that the strong coupling limit is in this 
regard a prototype of any Luttinger liquid, to recollect 
the lessons learned from the bosonic spin chain. There 
we learned that to break up squeezed space 'charge' fluc- 
tuations are needed changing its length from odd to even 
and vice versa. This implies that single charges can be 
created or annihilated and this is of course not a problem 
in a bosonic system because a single boson can condense. 
However, single fermions cannot condense and since in 
the Luttinger liquid for reasons just discussed the charge 
sector is fermionic, confinement is impossible. Admit- 
tedly, the argument is circular. It starts out postulating 
the existence of squeezed space as an entity unfrustrat- 
ing the spin system in the Marshall sign sense, to find 
out that the minus signs in turn offer a complete protec- 
tion of the squeezed space. This viewpoint suggests that 
there might be ways around the squeezed space and that 
states can be constructed which are beyond bosoniza- 
tion. Starting from strongly coupled microscopic dy- 
namics, one can image interactions which are strongly 
frustrating the spin system in the Marshall sign sense 
(i.e. longer range spin-spin interactions). Such interac- 
tions could lead to a 'signfuP spin physics in squeezed 
space, which in turn could diminish the 'statistical pro- 
tection', possibly leading to metallic states which are not 
Luttinger liquids. 

A final issue is, is there anything to be learned regard- 
ing the relevance of Luttinger liquid physics in higher 
dimensions? In this paper we have worked hard to per- 
suade the reader that squeezed space is a defining prop- 
erty of the Luttinger liquid. As such, it is a-priori not 
special to one dimension, in contrast to e.g. the lines of 
critical points and the Mandelstam construction. Given a 
complete freedom to choose the microscopic conditions, 
which fundamental requirements should be fulfilled to 
form squeezed spaces in higher dimensions? First, bi- 
partiteness is required and this is no longer automatic in 
higher dimensions. As a starting point one needs a Mott- 
insulator living on a bipartite lattice characterized by 
an unfrustrated, colinear antiferromagnet. Upon doping 



such a Mott-insulator the charges (holes) will frustrate 
this spin system unless special conditions are fulfilled: 
these holes have to form D — 1 dimensional connected 
manifolds as a fundamental requirement to end up in 
a bipartite space after the squeezing operation. Differ- 
ent from the one dimensional situation, true long range 
order will take over when it gets a chance. A first pos- 
sibility is that these D — 1 dimensional hole manifolds 
simply crystallize, forming charge ordered state accompa- 
nied by a spin system showing a strong ordering tendency 
as well, with the characteristic that the staggered order 
parameter flips every time a charge-manifold is crossed. 
One immediately recognizes the stripe phases which are 
experimentally observed in a variety of quasi-2D Mott- 
insulators, including the cupratesi Alternatively, assum- 
ing that the holes move in pairs, general reasons are avail- 
able demonstrating that the charge sector can turn into 
a superconductor (via a dual dislocation condensation 38 ) 
such that the manifolds continue to form domain walls 
in the sublattice parity although their locus in space is 
indeterminate. In direct analogy with the Haldane spin 
chain, such a state is characterized by an emergent 'sub- 
lattice parity' Zi gauge invariance. 

The above is just a short summary of some aspects 
of the 'stripe fractionalization' ideas and for a further 
discussion we refer to the literatura2*2i2iiS. Most im- 
portantly, the notion of squeezed space make it clear 
why 'Luttinger liquid-like' physics is not at all generic 
in higher dimensions but instead rather fragile, if it ex- 
ists at all. The bipartiteness of squeezed space-time in 
the space directions has to be protected and this requires 
microscopic fine-tuning. 

The punchline is that if one wants to contemplate man- 
ifestations of Luttinger liquid physics in higher dimen- 
sions it must be striped in one way or the other, since 
squeezed spaces are the most precise way to characterize 
the phenomenon of spin-charge separation as it arises in 
the specific one dimensional context. This insight also 
makes is clear why attempts to invoke the equations gov- 
erning the Luttinger liquids in whatever phenomenolog- 
ical spirit to explain physics in higher dimensions are 
bound to fail: these represent a dynamics which is slaved 
to an underlying geometrical principle which is only of 
the right kind in one space dimension. To bosonize the 
electron itself in two space dimensions one has to invoke 
geometrical/gauge principles of a fundamentally different 
kind 37 -22. 
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APPENDIX A: COMPUTATION OF THE 
CHARGE STRING CORRELATOR OF FREE 
SPINLESS FERMIONS. 



Ea. (|S5|) . We find that it can be fitted very accurately 
with the simple expression Eq. This may well be 

an exact result but we did not manage to find the solution 
with analytical means. 



In this appendix we discuss the numerical computa- 
tion of the free spinless fermion charge string operator, 



Using periodic boundary conditions, the charge string 
correlator can be written as 



((-1)£*=> C, ' ) ) = (k N ...h\(-l)^ 
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occupying the lowest N \k\ . . . k^) single fermion states. The product term on the last line equals —1 when yj S 

[l,x— 1] and 1 otherwise, taking into account the result of the factor (— 1)^j =1 spW . Part of this sum can be 
written as 
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abbreviating the second line with the 'star-delta function' 5*{p, k). Using this function, the expression (|A1() can be 
expressed as the determinant of a N x N matrix containing 6*(ki, kj) functions 



■SF(j)\ 



det 



/ 5*(ki,h) 5*(k 2 ,h) ■■■ S*(k N ,h) \ 

6*(h,k 2 ) S*(k 2 ,k 2 ) ••• S*(k N ,k 2 ) 

5*(Ai,fc3) S*(k 2 ,k 3 ) ••• S*(k N ,k 3 ) 

\s*(k u k N ) 6*{k 2l k N ) ••• 5*(k N ,k N )J 



(A3) 



and this determinant can be straightforwardly computed nu merically for a finite system. 



Careful analysis of the numerical data for a complete 
range of densities, demonstrates that, 




(A4) 

As an example, in Fig. we show results for 

Psf — y = ^ or A = 20 particles on a chain 
of length V — 200 and this compared with the 
analytic expression l)A4|l . In Fig. (|14|) the numeri- 

cal outcomes for the prefactor of ((— 1)^ =1 SFW } 



taking a normalization such that this prefactor is 
1 for psf = t? = 0-5- According to the exact re- 
sult by Parola and Sorella^ this prefactor is equal to 

( (_i ) E;: 1 1 »-o-) > ^z Bin( ^)/ ( ^^ c08(= ^ )) . 

The perfect match between this normalized numerical 
outcome and the function , 1 establishes the 

-y/sinOpsj?) 

density dependence of the amplitude in Eq. (|A4|I . 



In the thermodynamic limit V — > oo, y — > psf Eq. 
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(|A4|) becomes, 

( (-l)^ =1 )sf = 7== ■ 

(A5) 

reproducing the exact result by Parola and Sorella^ at 
the density psf = \- These authors showed that at this 
specific density the asymptotic form of D(x) is 

D(x) = ((-1)E;=o^« ) = A 2 V2 COSi ^Al \ (A6) 

\Jx + 1 
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FIG. 13: Numerical results (circles) for the function D{x) = 

((— 1)^-^=0 SF ^^ calculated from 1A3L as compared to the 
analytical form Eq. ljA4ft (full line) This is a representative 
example: we use N — 20 particles on on a chain of length 
V = 200 (density psf = S = 0.1), using periodic boundary 
conditions. 



APPENDIX B: A CAVEAT: SCHULZ' FLAWED 
LOGIC. 



density psf — — — 0.1. 
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FIG. 14: The crosses indicate the numerical results for the 
prefactor of the function D(x) = ((— 1)^-^=° SF ^') normal- 
ized to 1 for n = psf = f = 0.5 (V = 200). The full line 
corresponds with the function 1 

^/sin(vrp SF ) 

For completeness, let us list the outcomes for 

the expectation values (n(x)(— 1)^' =1 )s-F and 

( n(x)(— 1)^ =1 h(0)) which can be regarded as lat- 
tice derivatives of the charge string correlator Eq. i|A5() . 
Eg find 
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(A9) 



pendix we would like to comment on an argument due to 



(n(x) (-1)£J:> W ) SF 
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= sign cos(7rp S F - 1 y r—. = L - £ 5§ 1 , 

where the constant K is given by 

In addition, 

(n(x)(-l)^ nU) n(0)) SF = \[D{x-2)-2D{x-l) + D{ 

A 2 (co<&{tipsf) - 1) cos(irp SF 



x ) 
x) 



The idea that the charge fluctuates the space in which 
the spin system resides has a long history. In this ap- 
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the late Heinz Schula^S. His argument is not correct, but 
the flaw is subtle and informative regarding the workings 
of sublattice parity fluctuations. 

In the above, we re-derived the 'classic' result that the 
two point spin correlator (SS) ~ l/x K ° +Ks . Schulz 40 
asserted that this behavior can be explained by assum- 
ing that the system can be seen as a 1+1D harmonic 
crystal of charges in the continuum. The spins at the 
sites of this crystal would just form a Heisenberg anti- 
ferromagnet. True long range crystal order is impossible 
in 1+1D because the admixture of the Goldstone bosons 
(phonons) renders the correlations to be algebraic (alge- 
braic long range order, ALRO). Schulz' idea was simple: 
the spin systems does not live on fixed positions in space 



but instead on a medium undergoing gaussian fluctua- 
tions, as if the spin system 'surfs' on the gaussian charge 
waves. 

The effects on the spin correlator can be easily calcu- 
lated. In the continuum the spin density equals 

S( x ) = X! SHeis.{m)8(x - x m ), (Bl) 

m 

summing over all the electrons. Starting from the ALRO 
crystal, x rn can be written as x m = R m + u m , where 
R m = 222 ig the position in the m-th electron and u m its 
displacement. One finds for the correlation function, 



(S(x)S(0)) -Jd^ e-^^S Heis Xm)§ Heis Xm')}e^ R -- R ^^e^-^) . (B2) 

m,m' 

Due to the gaussian fluctuations, 



(e «?K,-« m ')) m \ m - m '\- a ^\ (B3) 



with a(q) ~ q 2 . The q integration in Eq. I)B2|) is dominated by the term q m irp = 2kp and using the Heisenberg 
correlation function, Eq. (|2T|l . 

(S(x)S(0)) ~ JdqJ2 e^ R -- R ^- x \S Hels Xrn)S Hels Xm'))\m^m'\- a( - 2k ^ 

p , ( i \m—m' 

« J <i<iY,<"'~ ~ 
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|m - m <y-+a(2k F ) 

(B4) 



This outcome looks indeed quite like the desired result, 
identifying a{2kp) with K c . However, this similarity is 
actually misleading! Schulz' crystal refers to the break- 
ing of translation symmetry by single electron charges. 
Implicitly, this refers to the strongly coupled regime con- 
sidered in the above and this crystal corresponds with the 
spinless-fermion ALRO crystal (e.g. reference filj) . The 
spinless fermion 2k p turns into a spinful electron 4k p 
wave vector. Accordingly, the exponent a(2kp) should 
be associated with the charge stiffness appearing in the 
4kp charge correlations, and this stiffness is not K c but 
instead 4K C \ For instance, in the large U case K c = 1/2 
and the Schulz argument would predict that the spin cor- 



relations would decay like 1/x 3 instead of 1/x 3 / 2 . 

Where is the flaw? In fact, the implicit assertion in 
the above is that (SS) ~ ( n n )ik F x (SS)Heis.- We 
learned, however, that the geometry of the spin system 
is fluctuated by kinks in their translational sector (the 
sublattice parity flips). These are dual to the charge- 
order and one has to use instead the exponentiated charge 

strings (SS) - (n(-l)£> n) x (SS)jieis.- As we showed, 
(-i)E« n) decays with an exponent which is K c itself. 
From the discussion in scction Hvl it is clear that this dual 
structure is in fact respected by bosonization. In this 
sense, bosonization 'knows' about squeezed space. 
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